Optical conductivity of the half-filled Hubbard chain 
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We combine well-controlled analytical and numerical methods to determine the optical conductiv- 
ity of the one- dimensional Mott-Hubbard insulator at zero temperature. A dynamical density-matrix 
renormalization group method provides the entire absorption spectrum for all but very small cou- 
pling strengths. In this limit we calculate the conductivity analytically using exact field-theoretical 
methods. Above the Lieb-Wu gap the conductivity exhibits a characteristic square-root increase. For 
small to moderate interactions, a sharp maximum occurs just above the gap. For larger interactions, 
another weak feature becomes visible around the middle of the absorption band. 
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In quasi one-dimensional materials like, e.g., organic 
semiconductors the electron-electron interaction ac- 
counts for the formation of excitons and, to a substantial 
degree, for the optical gap itself ||]. Therefore, a micro- 
scopic theory of the optical absorption in these materi- 
als requires a detailed understanding of one-dimensional 
Mott insulators 

The paradigm for a one-dimensional Mott insulator is 
the half-filled Hubbard model |] 
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It describes tt electrons with spin a =t:i which move 
from one site I of the carbon backbone to its neighboring 
sites. The lattice spacing is set to unity. The number of 
electrons N equals the number of lattice sites L. Using 
periodic boundary conditions, the kinetic energy is diag- 
onal in momentum space and gives rise to a cosine band, 
e{k) = —2tcos{k) of width W = 4t. The electrons' mu- 
tual Coulomb repulsion is mimicked by the purely local 
Hubbard interaction U . 

The model (|l|) is Bethe-Ansatz solvable The opti- 
cal gap is given by 
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For small and large U/t this gives A([/ < 2t) = 
{8t/TT)./U/icxp{-2nt/U) and A(;7 > it) = U - W + 
81n(2)t^/t/. The optical absorption is proportional to 
the real part of the optical conductivity which is related 
to the imaginary part of the current-current correlation 
fimction by ai{uj > 0) = Im{xjj(w > 0)}/uj, and 
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Here, |0) is the ground state, \n) are excited states, and 
Eq, En are their respective energies. Although -q — 0+ is 
infinitesimal, we may introduce a finite value to broaden 
our resonances at fiuj — En — Eq. In momentum space, 
J — ~{2et/K) X]/c (T sin(fc)c^o-^fc,<Ti is the current operator. 
We set Ti ~ \ throughout, and for our numerical results 
we use e = i = 1 in our figures. 

Eq. (||) shows why it is so difficult to calculate op- 
tical properties of Mott-Hubbard insulators. The spec- 
trum En is known exactly but very little is known about 



the oscillator strengths Tq „ — 



although it can 



be shown explicitly that they vanish unless En — Eq > 
A 0. Numerical calculations of (Ti(w) have been carried 
out using exact diagonalizations |^ and quantum Monte 
Carlo simulations Unfortunately, these approaches 
are seriously limited in accuracy or accessible system 
sizes. Since the system is a Mott insulator and not a 
Luttinger liquid metal, standard bosonization techniques 
cannot be applied to the half-filled Hubbard model JTo[ . 
Therefore, the calculation of the optical conductivity of 
this model is an important yet unsolved problem in the- 
oretical solid-state physics. 

In this work we employ the dynamical density-matrix 
renormalization group (DDMRG) method to de- 

termine the optical conductivity of the Hubbard insula- 
tor over the entire absorption spectrum. This numerical 
technique allows us to obtain ui (w) for all interaction 
strengths as long as the gap is not exponentially small, 
U > it. For large interaction strengths, U ^ t, we con- 
firm results obtained in the framework of a 1/U expan- 
sion jlSj. In the weak-coupling regime (A <S^t), we calcu- 
late ai{uj) analytically using exact field-theoretical meth- 
ods. The analytical results agree well with the DDMRG 
data for small U/t. 

We start our analysis with the large-?/ limit {U/t 
oo), where a rather simple band picture emerges p^ . If 
we ignore corrections of the order t/U electron transfers 
are limited to those processes which conserve the number 
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of double occupancies. Due to spin-charge separation, 
the oscillator strength Tg „ can be written as the convo- 
lution of the charge and spin contributions. The charge 
contribution follows from a simple band picture: we ex- 
cite one hole in the lower Hubbard band, eLHB(fc) = ^(k), 
and one double occupancy in the upper Hubbard band, 
euHB(fc) = U — e{k) (antiparallel bands). The total mo- 
mentum of the two charge excitations is q, and their en- 
ergy is oj. The spin sector enters the current-current cor- 
relation function via the momentum-dependent ground- 
state form- factor gq. 

For the large-?/ Hubbard model itself {U/t oo), the 
analysis of gq is rather involved. However, explicit analyt- 
ical results are available for the closely related problems 
of a dimerized (DIM) or a Neel-ordered (AF) spin ground 
state. For these states, gq vanishes unless g = or g = tt, 
g^^ = 2, ^ 0, and 5o°'^ = 9/4, g^^^ = 1/4 @. 
Now, optical transitions occur between two antiparallel 
bands {q = 0; eLHB(fc), euHB(fc)) and between two par- 
allel bands (g = tt; eLHB(fc)) euHB(^ + ti"))- Then, the 
optical conductivity becomes 
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For the large-?/ Hubbard model gq cannot be calcu- 
lated analytically. We can adopt a "no-recoil approxima- 
tion" to argue that the dominant contributions to 
the conductivity again come from q = and q = tt. Fig- 
ure |l| shows DDMRG results obtained on 64-site lattices 
for the large-C/ Hubbard model and the related models 
with a dimer or Neel-ordered spin ground state. On the 
scale of this figure there is no visible difference between 
our numerical results for the two models with spin order 
and the exact results with the same broadening r]. 
This shows that the DDMRG approach is very accurate 
in these limiting cases. 




FIG. 1. Current-current correlation function of the large-(7 
Hubbard model {U/t oo) and of two related models with a 
dimer and a Neel (AF) spin ground state {r] = 0.2). 



For the large-C/ Hubbard model DDMRG data are also 
in agreement with e g. (p| ) confirming the no-recoil ap- 
proximation in Ref. ]13|. In particular, one notes the 
small bump in ai{u;) at uj = U showing that g^^ > 
(see Fig. This bump is not visible in small chains but 
becomes more evident and sharper as we increase the sys- 
tem size. Therefore, there are two salient features in the 
optical absorption of the large-f/ Hubbard model, (i), the 
behavior at threshold is (Ti{uj) ~ y/uj ~ A, and (ii), the 
large density of states for excitations between parallel 
bands results in a small but visible peak in the middle of 
the absorption band. The second feature is not present in 
the optical conductivity calculated using a Neel-ordered 
spin ground state [|l3|,0 (see also Fig. |l]) because the 
spin form- factor gq of a Neel state vanishes for all g ^ 0. 

The density-matrix renormalization group (DMRG) is 
known to be a very accurate numerical method to de- 
termine static properties of low-dimensional lattice sys- 
tems Recently, Kiihner and White devel- 
oped an efficient scheme to calculate dynamical corre- 
lation functions in Heisenberg spin chains using DMRG 
and the correction vector method . We have extended 
this dynamical DMRG method to correlated electron sys- 
tems, such as the Hubbard model (|l|). Our method differs 
from that of Ref. Q in two points: (i) we calculate cor- 
rection vectors with Ramasesha's algorithm instead 
of a conjugate gradient method, and (ii) we compute dy- 
namical correlation functions directly from the correction 
vectors. 

Our DDMRG method allows us to calculate dynamical 
correlation functions, such as the r.h.s. of eq. (|3a|), very 
accurately for fairly large systems {L < 128) and a finite 
broadening factor rj. Thus, our numerical results always 
correspond to the actual function <7i{lo) convoluted with 
a Lorentzian L{llj) — ?7/[7r(aj^ + r/^)]. We compared the 
predictions of our method with exact results for a Peierls 
insulator |Q and the two models with a spin-ordered 
ground state (^) and found that errors were smaller than 
1% over the entire spectrum. In the limit of small U/t 
our numerical results agree very well with exact field- 
theory results, see below. We have also systematically 
checked various sum rules relating moments of the func- 
tion ai(ijj) to ground-state expectation values that are 
known exactly or can be evaluated with great accuracy 
using a ground-state DMRG method (l5[. For instance. 



doj 



t 

2L 



\J2{c+Ji+,^.+h.c.)\0) . (5) 



Lastly, the optical gaps deduced from DDMRG data al- 
ways agree with the exact result (|^) after taking finite- 
size corrections into account. 
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FIG. 2. Optical conductivity for U/t = 3, 6, 12 (from 
left to right), calculated with DDMRG on a 128-site lattice 
(77 = 0.1). Inset: ai{uj) for U/t = 12 (dashed) and 40 (solid) 
calculated on a 64-site chain {rj — 0.2). 

In Fig. 1^ we show the optical conductivity of the Hub- 
bard model calculated with DDMRG for several values 
of U/t. For U = AOt, (Ji{llj) resembles the large-?/ result, 
compare Fig. |l|. As the interaction strength decreases, 
the width of the absorption spectrum appears to dimin- 
ish slightly from 8t for the large-f/ Hubbard model to less 
than 7t for U — 3t. The small peak seen in the large-f/ 
Hubbard model is barely visible for U/t < 12, but this 
feature actually subsists at least down toU = At. As seen 
in Fig. ^, the shape of the optical conductivity changes 
progressively: for strong coupling {U — 12t), (Ji{ijj) dis- 
plays a broad distribution with a maximum clearly above 
the optical absorption threshold lo — IS.] for weaker cou- 
pling ([/ = 3<), (Ji(ijj) develops a very sharp peak close to 
the threshold and a long tail for higher frequencies. For 
U/t — !■ this peak eventually turns into a delta function 
as all the spectral weight is concentrated in the Drude 
peak for U — Q. 

For U/t > 3 the resolution of our calculations is good 
enough to see that the maximum of (7i(uj) occurs at 
a frequency significantly above the optical conductivity 
threshold. Our numerical data are compatible with a 
square-root increase of the optical conductivity above 
the Lieb-Wu gap for all values oi U/t > 3, in agree- 
ment with the analytical result for U ^ t. Obviously, 
the exact behavior of cri(w) at the threshold cannot be 
determined from the DDMRG results alone because the 
resolution is limited by the finite broadening rj and finite- 
size effects. For the same reasons, it is extremely difficult 
to study the optical conductivity of the Hubbard chain 
with DDMRG when the optical gap becomes exponen- 
tially small {U/t <3). 

For small interaction strengths, U <^ t, the half-filled 
Hubbard model can be mapped ||l^,|l^ onto its low- 
energy effective field theory, the SU(2) Thirring model 
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Here, i/j is a doublet of Dirac spinors and J° — ^tpj^a""!!) 
are SU(2) currents, where 7^ are 1-1-1 dimensional Dirac 



matrices and cr" are Pauli matrices. The theory (|^) de- 
couples into a massless (spin) and massive (charge) sec- 
tor. The charge sector exhibits a SU(2) symmetry as is 
required by the SO (4) symmetry of the half- filled Hub- 
bard model [ pO[ . This symmetry is in general broken 
in the Luther-Emery (a.k.a. U{1) Thirring) model pl| , 
and the optical absorption of the half-filled Hubbard 
model cannot be calculated using the exact solution at 
the Luther- Emery point |Q . 

In the field-theory limit, the current operator is found 
to be proportional to a SU(2) current j (x J dx Jf{x). 
This operator couples only to the gapped charge sector 
of the theory, so that only (multi) holon-antiholon scat- 
tering states contribute to (3b). It can be shown that lat- 
tice and band-curvature effects generate a coupling to the 
spin sector, but this is not important in the field-theory 
limit. We use the spectral representation of the excited 
states in terms of scattering states of holons (h), anti- 
holons (ft,) and spinous which form a basis of the Hilbert 
space P3| . Using charge-conjugation symmetry one finds 
that only states with equal number of holons Nh and 
antiholons Nj^ couple to . 

In the field-theory limit the holon/antiholon dispersion 
E{P) = \/P^ + (A/2)2 is parameterized in terms of a ra- 
pidity e as E{e) = (A/2)cosh6', P{e) = (A/2)sinh6'. 
A scattering state of N (anti)holons with rapidities 
{9j} and SU(2) indices {ej} {sj — h,h) is denoted by 
\9i,92, . . . ,dN)ei,e2.---eN ■ I^s energy and momentum are 
^^v = E{0j), Pn = Ef=i Pi^j)- In this basis, the 

two-point function of J fix) in the spectral representation 
reads 
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The matrix elements (0| (0)|6'i, . . . , On)ei...e_^ have been 
determined in ||2j], so that we can evaluate (Q). An inter- 
mediate state with Nh holons contributes to ui (w) only if 
Lo > NhA. Taking into account only intermediate states 
with Nfi — Nj^ = 1 we obtain (i^ = uj/A) 
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dx 1 — cos{x0 /it) cosh a; 
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where 9 — 2arccosh(j/). Formula (g) is exact in the in- 
terval A < u! < 2A. For u; > 2A there are corrections 
to (||), which arc due to multi holon/antiholon states and 
have a more complicated structure, but can be shown to 
be important only at energies uj ^ A ||2^. At present, 
the normalization C2(A) cannot be calculated analyti- 
cally. We note that S2p{cj) = cr2p(t^)/C2(A) is a univer- 
sal function of i' — lu / A. S'2p(A < a; < 2 A) is shown in 
the inset of Fig. |[ 
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FIG. 3. Optical conductivity for U/t = 3 from field theory 
(dashed) and DDMRG (solid). ai{ui) does not vanish below 
the gap A — 0.631t because a broadening of 77 = 0.1 is used. 
Inset: universal function S2p{^^) ~ (T2p(cij)/C2(A). 

(T2p('^) goes through a maximum for uj w 1.24A. 
For a; — !■ A we find the expected square-root behavior, 
(Tl(w A) ^ cr2p(^ ^ A) Vw" 
insulator fl^ , 



o'2p(w A) ~ y/uj — A. For the Peierls 
a square-root divergence occurs, reflecting 
the divergence in the density of states for excitations near 
the gap. For the Mott-Hubbard insulator, however, this 
divergence is suppressed by the momentum dependence 
of the oscillator strengths. For large U/t, eq. (^) shows 
explicitly that TQ.^{\k\ it, to ^ A) ^ sin^(fc) vanishes 
quadratically in the vicinity of the gap. 

To estimate the normalization C2(A), we use the ex- 
act sum rule (^. In the weak-coupling limit U <^t, most 
of the optical weight of the Hubbard model (0) must be 
concentrated at low energy w ~ A as (Ji{uj) reduces to 
a single Drude peak for C/ = 0. Therefore, in the field- 
theory limit, the l.h.s. of (||) can be determined from (||) 
as the omitted terms have a negligible contribution at low 
energy. For the r.h.s. we use the exact result at C/ = 0. 
We obtain C2(A) = 1.867 i/ A. 

In Fig. H we compare the field-theory prediction (0), 
using the exact value A « 0.63U of the optical gap ^) 
rather than its field theory value, and the DDMRG re- 
sult for U — 3t. We see that both results agree well up to 
Lu « 6A. This is surprising as the field theory is expected 
to work only as long as A <C It appears that even for 
U/t = 3, the contribution of states with more than one 
holon-antiholon pair to ai^Lu) is almost negligible in this 
range of frequencies. If we do the same comparison for 
larger U/ 1, we see that the range oiuj/A for which (T2p {^^) 
matches the DDMRG data becomes smaller, but both 
methods coincide around the optical absorption thresh- 
old even for fairly strong couplings {U — 6t). We have 
no explanation for this universal behavior of the optical 
conductivity above A. 

In conclusion, we have determined the optical conduc- 
tivity ai{uj) of the one-dimensional Mott-Hubbard insu- 
lator over the entire absorption spectrum for all values of 
the interaction strength. As also reported in we have 
found a square- root increase of 171(11;) at the absorption 
edge in contrast to the square-root divergence for the one- 



dimensional Peierls insulator. In principle, this difference 
could be used to distinguish experimentally between the 
two types of insulators. However, other interactions, e.g., 
the long-range parts of the Coulomb interaction, lattice 
structure, interchain couplings, and disorder, can signif- 
icantly change the absorption spectrum. We think that 
both dynamical-DMRG and field-theory approaches will 
permit the reliable calculation of dynamical properties 
for more general Hamiltonians which include these addi- 
tional interactions. 
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